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Abstract 

A leading order QCD analysis of spin asymmetries in polarized deep inelastic lepton 
nucleon scattering is presented within the framework of the radiative parton model. Two 
resulting sets of plausible leading order spin dependent parton distributions are presented, 
respecting the fundamental positivity constraints down to the low resolution scale = 
A^io ~ 0-23 GeV^. The dependence of the spin asymmetries is investigated 

in the range 1 < < 20 GeV^ and shown to be non-negligible for x-values relevant for 

the analysis of present data and possibly forthcoming data at HERA. 


The measurements of polarized deep inelastic lepton nucleon scattering yield direct 
information [1-3] on the spin-asymmetry 

A^ix 0^1 ~ ^ Q^) _ 

1 ^ ^ ^ ^ Ff^(x, Q2) Fi^(x, Q2)/ [2a;(l + R^(x, Q^))] ’ ^ ^ 

N = p,n and d = {p + n)/2, and R = Fl/2xFi = (F 2 — 2xFi)/2xFi, where subdominant 
contributions have, as usual, been neglected. The leading order (LO) QCD parton model 
relates A^{x,Q^) to the polarized {5q^,6q^) and unpolarized {q^,q^) quark distribu¬ 
tions, where q = u,d, s, via 


klf(a;,g^) = 


E.el 5q^{x,Q^)+6q^{x,Q^) 


Eg e2 [g^(a;, Q^) + q^{x, Q^)] 


( 2 ) 


due to the fact that 2gi = J2q 2Fi = J2q ^q{Q + q) and i? = 0 in LO (Callan- 

Gross relation). Henceforth we shall, as always, use the notation Sq^ = Sq and q^ = q. The 
experimental results are often presented and theoretically analyzed in terms of gi{x, Q"^) 
extracted via the measrtred values of A^^x, Q^), Ff{x, Q^) and R^{x, Q^) according to eq. 
(1). Frequently also the assumption about the independence of A^{x, Q"^) is employed 
in regions of x and where it has not been experimentally tested. This assumption is, 
as is well known Q, not theoretically warranted since one expects different Q^-evolutions 
of the numerator and the denominator due to the very different polarized and unpolarized 
splitting functions (except for 6P^g\x) = P^q\x)), respectively, especially in the small-a; 
region dominated by the flavor-singlet contributions. 


It should be further noted that an analysis of gi{x, as extracted from eq. (1) 
affords a full next-to-leading order (NLO) analysis due to the employed 7 ^ 0, 

typically of the order of 20 — 30%, and due to the fact that usually 2xFf{x, is well 
described only in a NLO QCD parton calculation since this latter quantity is typically 
about 20% smaller than the LO F:^. Such a NLO analysis is somewhat premature in 
view of the presently scarce and inaccurate data on gi {x, Q^). Some further ambiguities 
supporting this attitude are the ones concerning flavor SU (3) breaking effects, mentioned 
below, as well as those concerning the flavor singlet sector, in particular its gluonic com¬ 
ponent. For all these reasons it seems appropriate to apply hrst of all a LO analysis to 
the directly measured A^{x, Q"^) employing eq. ( 2 ), rather than to the derived gi {x, Q"^). 
Some further advantages of such an approach are that possible NLO and/or higher twist 
contributions are expected to partly cancel in the ratio of structure functions appear¬ 
ing in g^), in contrast to the situation for gf{x,Q‘^). Therefore we shall use all 
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presently available data [2,4-8] in the small-x region where ~ 1 GeV^ without bother¬ 
ing about lower cuts in usually introduced in order to avoid possible nonperturbative 
higher twist effects as mandatory for analyzing gi{x,Q‘^). It should be furthermore em¬ 
phasized that, within a consistent LO analysis {R^ = 0, i.e. 2xF^ = F.^), the results 
cannot be expected to agree with experimental extractions of g^{x,Q‘^) better than to 
about 20%. Alas, the extracted polarized LO parton distributions are also reliable only 
to within about 20%. The analysis affords some well established set of unpolarized LO 
parton distributions which will be adopted from ref. 0 . 


The searched for polarized LO parton distributions 6f{x, Q^), compatible with present 
data [2,4-8] on A^{x,Q‘^), are constrained by the positivity requirements implying 


\6f{x,Q^)\<f{x,Q^) 

where f = u,u, d, d, s, s, g, and furthermore by the sum rules 

Au + Au - Ad - Ad = gA = F + D = 1.2573 ± 0.0028 
Au +Au +Ad + Ad- 2{As + As) = 3F - D = 0.579 ± 0.025 

with the first {n = 1) moment Af defined by 

A/(q2 )^ f'dxdfix^Q^) 

Jo 


( 3 ) 


( 4 ) 

( 5 ) 


( 6 ) 


and the values of gA and 3F — D taken from |^. It should be noted that the first 
moment, i.e. the total polarization of each quark flavor Aq and Aq is conserved, i.e. Q^- 
independent as a consequence of helicity conservation at the quark-gluon vertex H, i.e. 
AP^q'^ = Jq dx6P^^\x) = 0 and AP® = 0. The constraint equations (4) and (5) are the 
ones used in most analyses performed so far. While the validity of the Bjprken ^fyi-sum-rule 
depends merely on the fundamental SU{2)f isospin rotation between matrix elements of 
charged and neutral axial currents, the constraint (5) depends critically on the assumed 
SU (3)/ flavor symmetry between hyperon decay matrix elements of the flavor changing 
charged weak axial currents and the neutral ones relevant for A/(Q^). We shall refer 
to the results based on the constraints (4) and (5) as the SU{3)f symmetric ‘standard’ 
scenario. 


There are some serious objections 0 to this latter full SU{3)f symmetry. As a 
plausible alternative, Lipkin [T^ has suggested a ‘valence’ scenario by assuming that 
the (flavor changing) hyperon /d-decay data fix only the total helicity of valence quarks 
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Aqv = Ag — Aq, i.e. 


gA = 1.2573 ± 0.0028 
3F-D = 0.579 ± 0.025 


( 4 ’) 

(S’) 


Ally — Ady — 

Auy + Ady = 

by assuming only the fundamental SU (2) f isospin symmetry u ^ d. We shall alternatively 
implement this extreme possibility in order to study the possible uncertainties, mentioned 
in the introduction, related to the SU{3)f symmetry breaking effects relevant for the 
^ n decay. 


To understand the important (practical) difference between the above two scenarios, 
let us consider the quantity 

f'dxgnx,Q^) ( 7 ) 

JO 

usually extracted [2,4-8] from measurements of Q^). In LO we have 

rr(Q^) = ±—Agg + —Agg +-AE 

1 vv ; 12 36 ^ 9 

= i^Agg + ^Agg + ^(As +As) (8) 

with the flavor nonsinglet (Agg^g) and singlet (AE) combinations being given by 


Agg = Au + Au — Ad — Ad 

Agg = Au + Au + Ad + Ad — 2(As -|- As) (9) 


AE = ^(Ag -|- Aq) = Agg -|- 3(As -|- As) 


In the ‘standard’ scenario Agg and Agg are entirely hxed by eqs. (4,5) which gives for Tg, 
for example, 

r; = l(F + z)) +A(3 f_£,) + 1(as + as) . (lo) 

For As = As = 0 we recover the original estimate of Gourdin and Ellis and Jaffe 0, 
^Iej - 0 .185. Therefore we need[| As < 0 in order to comply with recent experiments 
1^, H which typically give rf(Q^ = 3 GeV^) ~ 0.12 — 0.13. In the ‘valence’ scenario, only 
the valence contribution to Agg is hxed by eq. (5’), with the entire Agg still being hxed 
by (4’) due to the assumption Au = Ad = Ag, which gives for T^ in eq. (8) 


r? 


1(F + D) + ^{3F -D) + ^(lOAg + As + As-) 


( 11 ) 


^In a LO analysis the polarized quark first moments should be considered as effective moments since 
in NLO their detailed physical interpretation depends on the convention (factorization scheme) adopted 
where possibly [|^ Aq Aq — {as/‘iTr)Ag. 
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Thus, in contrast to eq. (10), a light polarized sea Aq < 0 will account for the reduction 
of rf as required by recent experiments [|, even for the extreme SU{3)f broken choice 
As = As = 0! We use this latter maximally SU{3)f broken polarized sea in our ’valence’ 
scenario, eqs. (4’) and (S’), which in addition is compatible with the SU{3)f broken 
unpolarized radiative input = 0 of ref. |^. A similar discussion holds for T” 

and of course also for Q"^). 


Apart from applying the above scenarios for the polarized input distributions to 
A^{x,Q‘^) rather than to gi{x,Q‘^), the main ingredient of our analysis is the imple¬ 
mentation of the constraint equation (3) down to [Q = 0.23 GeV^ which is 

not guaranteed in the usual studies done so far (recently, e.g. in [16-19]) restricted to 
A Qo ~ 1 ~ 4 GeV^. We follow here the radiative (dynamical) concept which resulted 
in the successful small-x predictions of unpolarized parton distributions as measured at 
HERA [^j. A further bonus of this analysis is the possibility to study the depen¬ 
dence of A^{x,Q‘^) in the small-x region over a wide range of which might be also 
relevant for forthcoming polarized experiments (HERMES) at HERA. The results for the 
Q^-dependence of Ai{x,Q^) are furthermore expected to hold almost unchanged also in 
the full NLO analysis due to the observed 0 perturbative stability of all the radiative 
model predictions for measurable quantities such as F^^x^Q"^). 


Turning to the determination of the polarized LO parton distributions 5/(x,(5^), it 
should be noted that 5g{x, Q^) does not appear explicitly in the LO expression (2) and 
is therefore only weakly constrained by present data [1,2,4-8] on A^{x,Q‘^) since it only 
enters indirectly via the Q^-evolution equations. It is thus necessary to consider some 
reasonable constraints concerning 6g{x, Q'^) in particular in the relevant small-x region 
as, for example, requirements of color coherence of gluon couplings at x ~ 0 (equal 
partition of the hadron’s momentum among its partons). This implies for the gluon and 


sea densities 20 




~ X 


as X 


( 12 ) 


where the scale Qo at which this relation is supposed to hold remains unspecihed. Al¬ 
though not strictly compelling, eq. (12) is expected [^ to hold at some ‘intrinsic’ bound- 
state-like scale (Qo ~ 1 GeV^, say), but certainly not at much larger purely perturbative 
scales 1 GeV^. Therefore we have htted our input distributions at Ql = = 0.23 

GeV^ using a general ansatz for the polarized (light) sea of the form ~ x"(l — x)^q 
with the result that all presently available asymmetry data require a ~ 0.9 — 1.1 for both 
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above scenarios. Consequently we have taken a = 1 and assumed this power also to hold 
for 6g{x, following eq. (12), although present data do not signihcantly constrain 

6g as will be discussed below. Our optimal LO distributions at = 0.23 GeV^ 

subject to the above constraints were found to be: 


6uv{x, ^i\o) 
5dv{x,g.\o) 
dq{x,g.lo) 
6s{x,filo) 


0.718x°'^My (a:, hIq) 
-0.728x°-^^dv{x, fxlo) 

—2.018a:(l — xY'^q{x, ii\q) 
6s{x,filo) = 0.72Sq{x,filo) 
16.55a:(l — xY'^‘^g{x, ii\q) 


(13) 


for the ‘standard’ scenario (corresponding to = 98.4/92 d.o.f.) respecting eqs. (4) and 
(5), while for the SU{3)f broken ‘valence’ scenario, based on the constraints (4’) and (S’), 
we have: 


Suv{x, p-lo) 
6dv{x,nlo) 

dq{x,gio) 

dg{x,g.lo) 


0.726x^'‘^^uv{x, pio) 
—0.668x^''^^dv{x, ii\o) 
—1.869a;(l — a;)°'^^g(a;, 

5s{x,g.\o) = 0 

Ux{l - xf-^^g{x, g^lo) 


(13’) 


which corresponds to = 96.8/92 d.o.f. The unpolarized input densities /(x, pIq) are 
taken from ref. and, for obvious reasons, we have not taken into account any SU{2)f 
breaking {Sit ^ 5d) as is apparent from our ansatz for Sq = 5u = Sd proportional to 
q = {it + d)/2 which should be considered as the reference light sea distribution for the 
positivity requirement (3). The fact that Ss{x, qt\o) 7 ^ 0 in (13) differs somewhat from our 
radiative input s(a:, pio) = 0; but for perturbatively relevant ~ 0.8 GeV^, where the 
leading twist-2 dominates in the small-a; region |P] , the positivity inequality (3) is already 
satished. In this respect the input (13’) for the ‘valence’ scenario, with the extreme SU{3)f 
breaking ansatz Ss{x,fi‘j^Q) = 0, is more agreeable as far as our radiative (dynamical) 
approach is concerned. Furthermore \Sqv{x, fJ>‘io)\ ~ qv{x,qt\o) as x ^ 1 in (13) and 
(13’) which is also compatible with arguments based on helicity retention properties of 
perturbative QGD [^. Finally, similarly agreeable hts to all present asymmetry data 
shown below (with a total of 101 to 103 for 92 data points) can be also obtained for 
a fully saturated (inequality (3)) gluon input Sg{x, = ( 7 (x, pIq) as well as for the 
less saturated 5g{x,gt\Q) = xg{x, qt^o). A purely dynamical input Sg{x, qt^^p) = 0 
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is also compatible with present data, but such a choice seems to be unlikely in view of 
Sq{x, f^io) 7^ 0) furthermore results in an unphysically steep being 

mainly concentrated in the very small-a: region x < 0.01, as in the corresponding case 
^ for the unpolarized parton distributions in disagreement with experiment. 


For calculating the evolutions of 5f{x, Q"^) to we have used the well known 


LO solutions in Mellin n-moment space (see, e.g, refs, p^), with the moments of the 
polarized LO splitting functions given in [|^. These were then Mellin-inverted to Bjprken- 
X space as described, for example, in |^. We have used / = 3 flavors when calculating 
and disregarded the marginal charm contribution to stemming from the sub¬ 
process 'y*g —>• cc [^. The LO running coupling utilized was a^lQ^) = 47r//3o ln((5^/A^) 
with ;3o = 11 - 2//3 and being given by |[ 


A(3.4,5) ^ 232, 200, 153 MeV . 


The number of active flavors / in as(Q^) was hxed by the number of quarks with 
taking rric = 1.5 GeV and rrih = 4.5 GeV. 


A comparison of our results with the data on Ai{x,Q‘^) is presented in hg.l. As 
already mentioned, £t results using a ‘saturated’ gluon 6g = g or 6g = xg at 


are very similar to the ones shown in hg.l. Note that {x, Q^) —*• const, as a; —»■ 1. The 
Q^-dependence of is presented in £g.2 for some typical hxed x values for 1 < 

< 20 GeV^. The predicted scale-violating Q^-dependence is substantial and similar 
for the two rather diherent input scenarios (13) and (13’). In the {x, Q^) region of present 
data [2,4-8], Af(x, Q^) increases with for x > 0.01. Therefore, since most present data 
in the small-x region correspond to small values of (~ 1 GeV^), the determination of 
gi{x, Q^) at a larger hxed (5 or 10 GeV^, say) by assuming Af(x, Q^) to be independent 
of Q^, as is commonly done [2,4-8], is misleading and might lead to an underestimate of gi 
by as much as about 20%, in particular in the small-a; region. The situation is opposite, 
although less pronounced, for —A^{x,Q‘^) shown in hg.2. This implies that \gi{x,Q‘^)\ 
might be overestimated at larger hxed by assuming ^"(a:, Q^), as measured at small 
Q^, to be independent of Q^. It should be emphasized that assuming Ai(a;, Q^) to be 
mdependent of contradicts, in general, perturbative QGD as soon as gluon and sea 
densities become relevant, due to the very diherent polarized and unpolarized splitting 
functions 0 , 6Pj;p{x) and Pip{x), respectively (except for 6P^^^ = P^^ which dominates 
only in the large-a; region). Moreover, the smaller x the stronger becomes the dependence 
of the exactly calculated Ai[x,Q‘^) on the precise form of the input at = /rio, as 
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can be seen in fig.2 for x = 10“^. For practical pnrposes, however, snch ambignities are 
irrelevant since, according to eq. (1), Ai(x,Q'^) ~ 2 xgi/F 2 —> 0 as x —> 0 is already 
nnmeasnrably small (of the order 10“^) for x ~ 10“^. Thns the small-x region is not 
accessible experimentally for gi{x,Q‘^), in contrast to the sitnation for the nnpolarized 
Fi 2 (x, Q^). It is interesting to note that the (approximate) asymptotic (x —> 0) QCD 
expression |^, Q for Ai(x, Q"^) does not even qnalitatively describe onr exact LO resnlts 
for Ai(x, in fig.2 for x > 10“^. 


In fig.3 we compare onr LO resnlts for gi {x^ Q^) with EMC, SMC and SLAC-E142/E143 
data as well as with the fit ’A’ of ref. [|^. Despite the fact that a LO analysis shonld 
not be expected to be more accnrate than abont 20%, as discnssed at the beginning, the 
EMC [|] and E143 p ‘data’ at fixed valnes of fall consistently below onr results in the 
small-x region. This is partly an artefact of the LO approximation (i.e. = 0 in eq. 

(1)) and partly due to the fact that the original small-x Af-data at small have been 

extrapolated ^ to a larger fixed value of by assuming Af(x, Q^) to be independent 
of Q^. According to the increase of Af with in fig.2, such an assumption underesti¬ 
mates gi in the small-x region at larger Q^. On the contrary, our results for do not 
show such a disagreement in the small-x region when compared with the SMC data 
in fig.3a where each data point corresponds to a different value of since no attempt has 
been made to extrapolate gi{x, Q"^) to a fixed from the originally measured A(^(x, Q^). 
Our predictions for the polarized LO parton distributions at the input scale 
in eqs. (13) and (13’) and at = 4 GeV^, as obtained from these inputs at 
for the two scenarios considered, are shown in figs.4a and 4b, respectively. The polarized 
input densities in fig.4a are compared with our reference nnpolarized LO dynamical input 
densities of ref. [Q which satisfy of course the positivity requirement (3) as is obvious from 
eqs. (13) and (13’). Our resulting polarized densities at = 4 GeV^ are compared with 
the ones (fit ‘A’) of ref. in fig.4b. Since the polarized LO gluon density 5g{x, Q"^) is 
not strongly constrained by present experiments, we compare our gluons at = 4 GeV^ 
in fig.5 with the ones which originate from imposing extreme inputs at Qq = such 
as 5g = g, 6g = xg and 6g = 0, instead of the one in (13’) for the ‘valence’ scenario. The 
results are very similar if these extreme gluon-inputs are taken for the ‘standard’ scenario 
in (13), and the variation of 6g{x,Q^) allowed by present experiments is indeed sizeable. 


Finally let us turn to the first moments (total polarizations) Af{Q‘^) of our polarized 
parton distributions, as defined in eq. (6). It should be recalled that, in contrast to 
Ag{Q‘^), the first moments of (anti)quark densities do not renormalize, i.e. are indepen- 
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dent of Q^, and thus the hrst moments implemented at the input scale = ^\q in (13) 
and (13’) remain the same at any Q^. Let us discuss the two scenarios in turn: 

‘standard’ scenario: From the input distributions in (13) one infers 

Auv = 0.9585, Mv = -0.2988, Ag = -0.0720, As = As = -0.0519, 
Agifilo) = 0.444, Ag(4GeV2) = 1.553, Ag(10GeV2) = 1.915, (14) 

which result in AS = 0.268. This gives, using eqs. (10) and (8), 

= 0.1506, F” = -0.0589 , (15) 

which, for a LO result, is in satisfactory agreement with recent SMG measurements 

[|,i 

F?(10GeV^) = 0.142 ±0.008 ±0.011, F"(5GeV2) =-0.08 ± 0.04 ± 0.04 (16) 

as well as with the most recent E143 data P implying F"(2GeV^) = —0.037 ± 
0.008 ±0.011. 

‘valence’ scenario: From the input distributions in (13’) one infers 

Auv = 0.9181, Ady = —0.3392, Aq = —0.0672, As = As = 0, 
Ag{iilo) = 0.417, Ag(4GeV2) = 1.509, Ag(lOGeV^) = 1.866, (14’) 

which result in a total singlet contribution of AS = 0.31. This gives, using eqs. 
(11) and (8), 

F? = 0.1478, F" = -0.0617 , (15’) 

which again compares well with the experimental values in (16). 

In both scenarios the Bjprken sum rule manifestly holds due to our constraints (4) 
and (4’), i.e. eq. (8) yields (A-u = Ad) 

r?(Q2) - ri(Q2) = ± . (17) 

6 

It is also interesting to observe that at our low input scale = h\q = 0.23 GeV^ the 
nucleon’s spin is, within 20%, carried by quarks and gluons, |AS ± Ag^gi^o) k. 0.57, 
according to (14) and (14’), which implies for the helicity sum rule 


1 = IaE + A9(Q^) + L,XQ^) 


(18) 


Lz{fi\o) ~ 0. The approximate vanishing of this latter nonperturbative angular momen¬ 
tum, being built up from the intrinsic kx carried by partons, is intuitively expected for 
low (bound-state-like) scales but not for ^ ^\q. From our results for the total gluon 
polarization Ag{Q‘^) in (14) and (14’) it is also apparent that the Q^-independent anomaly 
contribution —[aaiQ"^)/‘^T^]Ag{Q‘^) ^ —0.03 could equally well serve^ as the exper¬ 
imentally required negative contribution to F)^ in (8), instead of allowing for hnite 
negative sea contributions As and Ag in eqs. (10) and (11), respectively. Such alternative 
factorization scheme scenarios cannot be decided on purely theoretical grounds, but nev¬ 
ertheless a consistent and convention independent analysis of the anomaly contribution in 
Bj0rken-a; space requires the full knowledge of all polarized two-loop splitting functions. 

Finally let us remark that our results demonstrate the compatibility of our very re¬ 
strictive radiative model, cf. eq. (3), down to = pIq = 0.23 GeV^, with present 
measurements of deep inelastic spin asymmetries. A Fortran package containing our 
optimally htted ‘standard’ and ‘valence’ distributions can be obtained by electronic mail 
from vogelsang@v2.rl.ac.uk . 
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Note added 


While this manuscript was being completed, a complete calculation of the two-loop split- 

,{i) 


ting functions {x) in the MS factorization scheme has appeared for the hrst time 


For the reasons stated in the introduction we postpone a full quantitative NLO analysis 
to a future separate publication. 
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Figure Captions 


Fig.l Comparison of our results for Ai{x,Q‘^) as obtained from the fitted inputs at 
= h\q for the ’standard’ (eq.(13)) and ’valence’ (eq.(13’)) scenarios with present 
data [2,4-8]. Our results at different values of x correspond to different values of 
according to experiment where each data point refers to a different value of 
starting at ~ 1 GeV^ at the lowest available x-bin. 

Fig.2 The dependence of Q"^) as predicted by the LO QCD evolution at various 

hxed values of x. 


Fig.3a Comparison of our ’standard’ and ’valence’ scenario results with the data [2,4-6] 
for gl’’^{x^Q‘^). The SMC data correspond to different ~ 1 GeV^ for x > 0.005, 
as do the theoretical results; here the comparison with the GS £t [^ is limited to 
X > 0.01 since the data at lower x correspond to < 4 GeV^. The GS £t result 
for the E143 {gi) data corresponds to a hxed = 4 GeV^. 

Fig.3b Same as in fig.3a but for g'^{x^Q‘^). The E142 and E143 data [^ correspond 
to an average (Q^) = 2 and 3 GeV^, respectively. 

Fig.4a Comparison of our htted ’standard’ and ’valence’ input densities in eqs. (13) and 
(13’) with the unpolarized dynamical input densities of ref. P]. 

Fig.4b The polarized densities at = 4 GeV^, as obtained from the input densities 
at = jji^Q in £g.4a. The htted GS (set A) densities of ref. are shown for 
comparison. 


Fig.5 The experimentally allowed range of polarized gluon densities at = 4 GeV^ for 
the ’valence’ scenario with diherently chosen 5g{x, inputs. The ’htted 5g'’ curve 
is identical to the one in hg.4b and corresponds to dg^x^n^o) in eq. (13’). Very 
similar results are obtained if 5g{x,gi\Q) is varied accordingly within the ’standard’ 
scenario. 
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Fig. 3b 
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